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ABSTRACT 

A quantitative study of the Luminous Blue Variable NGC 2363-V1 in the Magellanic galaxy 
NGC 2366 (D = 3.44 Mpc) is presented, based on ultraviolet and optical Hubble Space Telescope 
STIS spectroscopy. Contemporary WFPC2 and William Herschel Telescope imaging reveals a 
modest V-band brightness increase of ~0.2 mag per year between 1996 January-1997 November, 
reaching V=17.4 mag, corresponding to My—— 10.4 mag. Subsequently, VI underwent a similar 
decrease in V-band brightness, together with a UV brightening of 0.35 mag from 1997 November 
to 1999 November. 

The optical spectrum of VI is dominated by H emission lines, with Fen, Hei and Nai also 
detected. In the ultraviolet, a forest of Fe absorption features and numerous absorption lines 
typical of mid-B supergiants (such as Sill, Si in, Siiv, Cm, Civ) are observed. From a spec- 
tral analysis with the non-LTE, line-blanketed code of Hillier & Miller (1998), we derive stellar 
parameters of T,=llkK, i?»=42Oi? , log (£/L )=6.35 during 1997 November, and T*=13kK, 
i?*=315i?Q, log (L/Lq)=6A for 1999 July. The wind properties of VI are also exceptional, with 
M ~ 4.4 x 10~ 4 M Q yr _1 and ~300 kms -1 , allowing for a clumped wind (filling factor = 
0.3), and assuming H/He~4 by number. 

The presence of Fe lines in the UV and optical spectrum of VI permits an estimate of the 
heavy elemental abundance of NGC 2363 from our spectral synthesis. Although some deficien- 
cies remain, allowance for charge exchange reactions in our calculations supports a SMC-like 
metallicity, that has previously been determined for NGC 2363 from nebular oxygen diagnostics. 

Considering a variety of possible progenitor stars, VI has definitely undergone a giant erup- 
tion, with a substantial increase in stellar luminosity, radius, and almost certainly mass- loss rate, 
such that its stellar radius increased at an average rate of ~4 kms during 1992 October - 1995 
February. The stellar properties of VI are compared to other LBVs, including 77 Car and HD 5980 
during its brief eruption in 1994 September, the latter newly analyzed here. The mass-loss rate 
of the HD 5980 eruptor compares closely with VI, but its bolometric luminosity was a factor ~6 
times larger. 

Subject headings: stars: early-type - stars:mass-loss - stars: variables: other - stars: individual: 
NGC 2363-V1 



1 Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, 
which is operated by AURA, Inc., under NASA contract NAS5-26555. 
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1. INTRODUCTION 

Luminous Blue variables (LBVs) are evolved massive stars which display photometric and spectroscopic 
variability with different timescales and amplitudes. Although LBVs show micro-variations (Amag ~ 0.1- 
0.3 mag; At ~ days-months) like other supcrgiants, it is their moderate "excursions" (Amag ~ 1-2 mag at 
constant M),oi; At ~ years) and giant eruptions (Amag > 3 mag with an increase in M),oi; At ~ years - 
decades) that define them as a class. They are found in the upper left section of the H-R diagram, having 
Mboi < —9 and log (T e //) > 4.1 (except during an eruption when their surface temperature often drops 
to log (T e ff) — 3.85). LBVs are often, but not always, surrounded by a small (1-3 pc) nebula resulting 
from a giant eruption. We refer the reader to Humphreys & Davidson (1994) and Nota & Lamers (1997) for 
comprehensive, recent reviews on the subject. Although it is now widely accepted that the LBV phenomenon 
represents a normal yet unstable evolutionary phase in the life of the most massive stars, the exact mechanism 
leading to the major eruptions is still unknown. This lack of understanding is due in part to the problems 
involved in the theoretical modeling of the eruption phenomenon (Stothers 1999b), but also to the rarity 
of such events which has until recently prevented us from determining the physical parameters of an LBV 
during an eruption. 

In 1996, we reported the discovery of a bright variable star in the giant extragalactic H n region 
NGC2363 (Drissen, Roy & Robert 1996, 1997, hereafter DRR97). This star, known as NGC2363-V1 (here- 
after VI), was first noticed in HST/WFPC2 images. An archival search for ground-based images allowed us 
to determine that this star was below the detection limit in 1991 and 1992 (V > 21.5), but became visible 
in late 1993. In early 1995, it became the brightest star in its galaxy, with V — 18.0 mag. The photometric 
behavior, absolute magnitude and the presence in its spectrum of a strong Ha emission line suggested that 
VI is an LBV currently experiencing a major eruption event. 

The host of the current LBV event is the giant Hn region NGC 2363, located at the end of the bar in 
the Magellanic galaxy NGC 2366 (D = 3.44 Mpc; Tolstoy et al. 1995). The mctallicity of NGC 2363 has 
generally been determined to be close to that of the SMC (Gonzalez-Delgado et al. 1994), but Luridiana, 
Peimbert & Lcitherer (1999) suggest that the metallicity of this region has been underestimated, and that 
a value of Z = 0.25Z Q would give a better fit to the observed line ratios. A detailed study of the stellar 
content of NGC 2363 (Drissen et al. 2000) reveals that two clusters of different ages are responsible for the 
ionization of the nebula. The main contributor is a very young (1 Myr or less), dense cluster of massive stars 
still embedded in their natal molecular cloud. VI is a member of the second cluster, which is between 3 and 
5 Myr old and contains 3 Wolf-Rayet stars. 

Given the importance of this star in our understanding of the LBV phenomenon, we proposed a spec- 
troscopic and photometric follow-up of VI with the Hubble Space Telescope. This project was granted the 
long-term status for cycles 7 to 9 (GO- 7391, 8403 and 8781). We present in this paper an analysis of 
HST/STIS spectrograms obtained in 1997 November and 1999 July, covering the wavelength range 1180 A 
to l,um, supplemented by ground-based and HST photometry. We present our new observations in § 2 and 
produce a quantitative analysis with the non-LTE line blanketed code of Hillicr & Miller (1998) in § 3. A 
plethora of Fe n lines in the UV and optical spectrum of VI permits an estimate of the heavy metal content 
for NGC 2363. Our results are discussed in § 4, with conclusions reached in § 5. 
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2. OBSERVATIONS 

DRR97 published a light curve of NGC 2363-V1 between 1991-1996 revealing a dramatic increase of 
>3.5 mag in brightness, although no color information was available. In this section, we present recent 
HST/WFPC2 and ground-based imaging, plus HST/STIS spectroscopy. 

2.1. Photometry 

UBVRI images of VI were obtained on the night of 1998 January 14/15 with the 4.2 m William Herschel 
Telescope (WHT) on La Palma, Canary Islands. The Cassegrain auxiliary port was used with a Tektronix 
1024 x 1024 CCD (O.'Tl/pixel). The exposure times were 300 s at BVRI and 600 s at U. The seeing was 
measured to be 0.9-l".0, sufficient to resolve VI from the nearby Hn region. Landolt (1992) standard star 
fields (98 624, 98 626 and 98 634) were also observed to provide the zero point calibration for the VI images. 

Accurate magnitudes for VI are impossible to obtain using conventional aperture photometry because of 
severe contamination from the nearby H n region. The following approach was therefore adopted. Magnitudes 
were measured using a 35 pixel aperture with DAOPHOT (Stetson 1987) for several isolated stars in the 
same frames as VI. The same stars and VI were then measured using a much smaller aperture of 3.5 
pixels, and sky inner and outer radii of 4 and 7 pixels. These values were chosen to minimize the Hn region 
contamination. The UBVRI magnitudes of VI were then determined by assuming that the same fraction of 
light has been excluded from the smaller apertures for all the measured stars. 

We also obtained WFPC2 photometry with the F170W, F336W, F547M and F1042M filters on 1998 
March and December, 1999 November and 2000 February. The images were processed with the standard 
pipeline and aperture photometry of VI was performed with DAOPHOT. The prescriptions given by Whit- 
more (1995) and Holtzman et al. (1995) were then followed to transform the instrumental magnitudes into 
the Johnson system. 

Finally, near-infrared imagery of NGC 2363 were obtained in 1997 January by D. Devost and R. Doyon 
with the Monica Camera attached to the Canada-France-Hawaii telescope, which is described in more details 
in Drisscn et al. (2000). From these observations, we get J=17.6 (±0.1), H = 17.2 (±0.1) and K = 16.9 
(±0.2). These values are slightly fainter than the predicted magnitudes based on the models described in 
the next sections (J = 17.1, H = 16.9 and K = 16.7), indicating that cold dust has not yet formed in the 
stellar outflow. 

The results of all photometric measurements are given in Table 1. The agreement between the ground- 
based and HST March data (taken 6 weeks apart) is very good in the V band, but the difference in U (0.6 
mag) is significant. An updated light curve with the new data (in the V band) is shown in Figure 1. After 
the original outburst in 1993-1995, the brightness increase has been much more modest, though constant 
during 1996-1997, with a rate of ~ 0.2 mag per year, as if the star had reached a state of quasi-equilibrium. 
The 1998-2000 photometric observations suggest a small, but comparable, faintening in the V band. The 
WFPC2 data also indicate that VI brightened by ~ 0.3 mag at 170 nm between 1998 March and 2000 
February; during the same period, it faded by about the same amount at 1 /um. 
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2.2. Spectroscopy 

We obtained HST/STIS spectra of NGC 2363-V1 in 1997 November and 1999 July, using the MAMA 
detectors short ward of 3000A and the CCDs at longer wavelengths. Because of the sensitivity of the MAMA 
detectors to high radiation levels in the South Atlantic Anomaly, the 1997 optical and UV spectrograms of 
VI had to be acquired on separate occasions three weeks apart (November 2 and 21). Given the relative 
photometric stability of VI between 1997 November and 1998 December, we do not expect that the spectro- 
scopic appearance of the star has significantly changed during this short interval. In 1999, the optical and 
UV spectra were obtained 2 days apart. The following gratings were used: G140L (160 minutes) and G230L 
(90 minutes) with the MAMA detector; G430L (50 minute exposure), G750L (35 minutes) and G750M (Ha, 
medium resolution; 50 minutes) with the CCD detector. The 0.2 arcsecond-wide slit was used in all cases, 
which corresponds to a linear scale of 3 pc at the distance of NGC 2366. 

Because the data originally processed with the standard HST pipeline showed unexplained anomalies, 
we have retrieved our data from the Canadian Astronomical Data Center (CADC). The CADC archives uses 
the best available calibration files to re-calibrate the data. Moreover, because VI is located in the NGC 2363 
nebula, extra care has been taken in subtracting the adjacent nebular spectrum from the spectrum of VI. 
Nebular contamination is small (VI is in the middle of an expanding cavity) and relatively easy to subtract 
in the high (spatial) resolution STIS data compared with the highly contaminated ground-based spectra, 
but the standard pipeline calibration did not do a good job in this respect. Standard IRAF procedures 
were used. Despite the use of contemporary flatfield calibrations, the far red portion of the spectra is 
contaminated by fringing. By convolving our observed spectrophotometry with suitable synthetic filters, we 
obtained wide-band Johnson photometry (see Table 1) that is in excellent agreement with WFPC2 imaging. 

The low- resolution flux calibrated HST datasets are presented in Figure 2. Balmer emission lines are 
the most prominent spectral features in the visible; weaker lines at Hei A4471, 5876, 6678, Nai A5890 plus 
numerous Fen transitions are also observed. The apparently noisy near-UV depression results from many 
Fe-transitions, rather than high interstellar reddening. 

Since individual features are generally unresolved with the low resolution gratings, we have also obtained 
higher resolution Ha STIS observation with the G750M grating for each epoch. These arc presented in Figure 
3 and reveal a very strong Ha P Cygni profile betraying the extent of VI 's envelope and its high mass-loss 
rate. The observed Ha flux is 5.6xl0~ 14 and 6.3xl0~ 14 crgs^cm- 2 for 1997 November and 1999 July, 
respectively. Correcting for interstellar reddening (see below) and the distance to NGC 2363, the intrinsic 
Ha luminosity of VI is 2.4-2. 7x 10 4 L Q . 

While the general appearance of the spectrum is similar at both epochs, a few changes are worth noting: 

(1) The UV continuum flux (1200 - 2800 A) is stronger in 1999. However, the line fluxes in this 
wavelength range has not changed significantly. 

(2) The strength and profile of the Balmer lines has changed (see Table 2). The flux of the Ha line in 
the July 99 dataset is ^20% stronger and its blue absorption component extends to higher velocities (Figure 
3). On the other hand, the equivalent widths of the other Balmer lines have slightly decreased and the 
absorption components of their P Cygni profile are less deep in 1999. 

(3) The Fell lines at 4925 and 5016 A which were strong in the 1997 data became much weaker in 1999. 
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2.3. Morphology of the UV spectrum 

The resolution of the STIS data (~ 300 km s _1 ) does not allow us to easily, and unequivocally, distinguish 
between the interstellar and purely stellar contribution components to the resonance lines on the basis of 
their width alone. But a direct comparison of the line strengths between the spectrum of VI and those of 
knots NGC 2363-A and B (obtained with HST/FOS; see Drissen et al. 2000 and Figure 4 here) should allow 
us to identify the lines which are intrinsic to the star's photosphere. We assume that since VI and knot B arc 
both located within the same expanding supcrbubble, the spectral contribution from their local ISM will be 
similar, although knot B is a cluster of OB stars and some ISM lines in its spectrum could be contaminated 
by a stellar contribution. On the other hand, the ISM local to knot A could be different, but contributions 
from OB stars are negligible. The contribution from Milky Way gas should be the same in all three objects. 
From Figure 4, which shows the rectified spectrum of VI in 1999 in the range 1160 - 1700 A along with that 
of knots A and B, and from Table 3, which lists the equivalent width of the most conspicuous UV lines in the 
three objects, we conclude that the UV lines of VI have a strong photospheric contribution. As examples, 
one may compare the lines of Sin A1260, Cn A1335, Siiv A1400 and Sin A1527. None of these UV lines in 
VI however show any evidence for an emission component, suggesting a rather low effective temperature. 

The following lines are both strong in VI and absent or very weak in the spectrum of knots A and B; 
they arc thus likely to have a photospheric origin: Cm A1175, Mgn A1240, Sn A1251,54,60, Sin A1265, 
Tim+Silll between 1286 and 1303A, Sin A1309, and Sin A1533. 

The Cm A1175 line typically shows a strong P Cygni profile in late- type O to early- type B supergiants 
and becomes mostly photospheric in mid B supergiants (Snow & Morton 1976). It is an obvious iW\ = 2.4 
A) absorption feature in VI, but appears much weaker than in late-type B supergiants such as rj CMa (B5Ia) 
or (3 Ori A (B8Ia). Also, features seen in the spectrum of VI between 1295 and 1303 A, are often associated to 
Sim in B stars (Massa 1989). Therefore in Figure 4, we have plotted along with VI the rectified spectrum of 
an average Galactic B5.5I star (based on high resolution IUE spectra of 4 B6I and B5I stars from the library 
of Robert 1999). The resolution of this B5.5I spectrum was degraded to mimic the STIS data. Similarities 
with VI are striking. In the following discussion we compare the spectral feature of VI with those of B stars. 

The strength of the Si n A1265 line is a strong (inverse) function of the effective temperature in OB stars, 
but insensitive to the luminosity class (Massa 1989; Prinja 1990), increasing from 100 mA at T e ff — 25000K 
to 1600 mA at T eff = 13000K. Sin A1265 is not seen in stars hotter than B0. The large value of W\ (Sin 
A1265) = 1800 mA in the 1999 spectrum of VI sets a relatively stringent upper limit of T e ff < 13000K, 
albeit based on a Galactic supcrgiant calibration. 

The strong line at ~ 1303 A is a blend of stellar and interstellar lines. The interstellar component 
(mostly Oi A 1302 and Sin A1304) dominates in knots A and B (total W\ = 2.2 A in both cases), but the 
blend is stronger in VI (W\ = 4 A), indicating a strong photospheric contribution from Sim A 1301,03 and 
Sin A1304 typical of late to mid-B stars. 

The Sin A1309 line behaves like Sin A1265 as a function of OB spectral types. The Sin A1533 is more 
difficult to study in hot OB stars as the broad wind profiles developing in Civ A1550 engulfs it and as it is 
blended with Feiv (Nemry et al. 1991) in hotter stars. Nevertheless there are indications that Sin A1533 is 
also a good B star signature. 

Spectral features around 1317, 1345, and 1455 A, which are associated with Nin, Pm, and Tim, are 
typical of B supergiants; they are not seen in giants and dwarfs B (except maybe in the earliest types B0-1). 
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The strength of the doublet Si iv A1394,1402 is strongly temperature and luminosity dependent (Walborn 
et al. 1995). It is absent from all late B luminosity classes, appears in absorption in mid-early types, and 
in the case of the early B supergiants shows a strong P Cygni wind profile. An absorption doublet of C iv 
A1549,51 is only seen in hot BO-2 dwarfs and giants. A wind profile appears in the hottest giant BO. In the 
case of supergiants, the absoption is seen for mid B and quickly becomes a strong P Cygni profile in type B4 
and hotter. We do not see Si iv and C iv wind profiles in VI but rather absorption typical of B5-6 I stars. 

Around 1600 A we see many strong lines of Fell and Fein which resemble the features of B supergiants. 
These lines are much weaker in B giants and dwarfs. 

We therefore conclude that the UV spectrum of VI is similar to that of a B supergiant star of type 5-6. 

3. SPECTROSCOPIC ANALYSIS 

An important consideration for the interpretation of Vl's spectral energy distribution is the determi- 
nation of the extinction. Gonzalez-Delgado et al. (1994) derived E(B-V)~0.15 mag from a ground-based 
integrated nebular spectrum of the whole H n region. Since the extinction is likely to be different in knot 
A, which dominates the global nebular spectrum, than in the middle of the expanding bubble where VI and 
knot B lie, we have extracted the spectrum of the nebula from the pixels adjacent to VI in the STIS long- 
slit spectroscopy. The nebular Ha/H/3 line strengths from this region reveal a low value of E(B-V)=0.06 
using Case B recombination theory (Hummer & Storey 1995); this is in perfect agreement with Drissen 
et al. (2000) who recently estimated E(B-V)=0.04 (Galaxy) + 0.02 (NGC 2363) based on HST/FOC UV 
spectra of knot B in NGC 2363. Therefore, if an associated circumstellar nebula is present, it does not add 
measureably to the interstellar extinction, arguing against a dense, rj Car- like nebula surrounding VI. We 
have thus reddened the theoretical continuum energy distributions by E(B-V)=0.04 using the Scaton (1979) 
Galactic extinction law, plus E(B-V)=0.02 using the Bouchet et al. (1985) SMC law. Assuming a distance 
of 3.44Mpc to NGC 2366 implies M v =-10.44 mag during 1997 November, and -10.29 in 1999 July. 

For stars with extended atmospheres such as NGC2363-V1, as evidenced from its P Cygni Ha profile, 
the usual assumptions of plane parallel geometry and local thermodynamic equilibrium (LTE) are totally 
inadequate. In this section we introduce the atmospheric model that we shall use for VI and discuss our 
analysis technique. 

3.1. Atmosphere code 

We utilize the iterative, line blanketed technique of Hillier & Miller (1998) which solves the transfer 
equation in the co-moving frame subject to statistical and radiative equilibrium in an expanding, steady-state 
atmosphere. Populations and ionization structure are consistent with the radiation field. Line blanketing is 
treated via a global Doppler line width of Vbop^lO kms" 1 . 

Despite the high quality dataset available, the spectroscopic analysis of NGC 2363-V1 has proved to be 
difficult, principally because its low excitation state means that Fen is the dominant blanketing ion, which 
is extremely complex. In addition, as discussed by Hillier et al. (1998), mass-loss and hydrogen content are 
highly coupled. Since we are unable to tightly constrain the stellar temperature of VI, due to the absence 
of lines from adjacent ionization stages of the same element, we cannot derive a mass-loss rate without first 
fixing its He abundance. Consequently, we shall adopt an abundance of H/He=4 by number, which is typical 
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of LBVs (Crowther 1997) , and subsequently investigate the effect of different He abundances on the mass-loss 
rate. Similarly, metal abundances were fixed at 0.20Z Q as indicated by Gonzalez-Delgado et al. (1994) for 
NGC 2363, and varied once stellar parameters were derived. 

A simplifying 'super level' approach is used for individual levels (Anderson 1989), particularly for iron- 
group elements. In this approach, several levels of similar energies and properties are treated as a single 
'super level', with only the populations of the super level included in the solution of the rate equations. The 
population of an individual atomic level in the full model atom is determined by assuming that it has the 
same departure coefficient as the corresponding super level to which it belongs. 

The present model include representative model atoms for H, He, N, Mg, Al, Na, Ca, Si and Fe, as 
indicated in Table 4, in which a total of 32,191 transitions are considered. Details of each ion are included, 
such that for hydrogen, n=l to 30 full levels are considered, which are grouped into 10 super levels. Our 
default atomic dataset for Fen was taken from Nahar (1995), in which 827 atomic levels, were combined 
into 134 super levels, producing 21,907 transitions (22,924 weak transitions with log gf < —4 were omitted). 
Hillicr et al. (2000) have recently analyzed the UV and optical spectrum of i] Car based on the Nahar 
compilation and found improved comparisons with observation if Kurucz (1993) oscillator strengths were 
selected instead of those from Nahar, although negligible differences were obtained for VI. 

In the outer wind of VI, and other LBVs, hydrogen becomes neutral, so it is necessary to allow for charge 
exchange reactions. Of particular relevance to LBVs is the charge exchange reaction: Fe 2+ 3d 6 ( 5 D) + H 
<-> Fe + 3d 6 4s (a 6 D) + H + (Neufeld & Dalgarno 1987) since the ratio of charge recombination to radiative 
recombination for Fe 2+ is ~ 10 3 H/H+ (Hillier et al. 2000). This has a major influence on the ionization 
structure, and significantly enhances the Fe n line strengths. Also relevant to LBVs is the nitrogen reaction: 
N 2 + + Hc+ <-> N+ + Hc 2 + (Hcrrcro et al. 1995). 

Finally, since it is well established that the dense winds of early-type stars are clumped (e.g. Lepine & 
Moffat 1999; Eversberg et al. 1998), we consider a relatively simple volume filling factor approach for VI 
following Hillier & Miller (1998). Estimates of clumping factors in the optical line forming region for VI are 
provided by simultaneously matching recombination emission line strengths (sensitive to the square of the 
density) and electron scattering wings (sensitive to density) as discussed by Hillier (1991). 

3.2. Stellar parameters 

Since we were unable to use lines from adjacent ionization stages of one particular element to derive 
stellar temperatures, we instead relied on the stellar continuum, de-reddened as described above. We have 
calculated a grid of models for VI at each epoch in order to simultaneously match the stellar continuum 
distribution from which the stellar temperature and luminosity are derived, plus Ha which provides an 
excellent tracer of mass-loss in the winds of early-type stars (Puis et al. 1996). We have adopted a typical 
H/He content of LBVs for NGC2363-V1, namely H/He=4 by number (see also below). 

For the 1997 November dataset, we find optimum agreement for a model in which T» = llkK, log(L/L0)=6.35 
and _R(tr oss = 10) = 420i? Q , as shown in Figure 5 (upper panels). Since the atmosphere of NGC2363-V1 
is very extended, -R(tr oss = 2/3) = 590i? Q , the formal 'effective temperature' is T e ff=9,200K. Similar prop- 
erties for July 1999 are T» = 13kK, log(L/L Q )=6.4 and R(tr oss = 10) = 315i? Q , shown in Figure 5 (lower 
panels). Therefore, the spectroscopic differences between the two epochs represent variations in stellar radius 
(i.e. temperature) at fairly constant bolometric luminosity. This naturally explains the simultaneous UV 
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brightening and optical fading during this period. 

While the agreement between the model and the observations is generally very good, we wish to highlight 
several important discrepancies revealed in this work. The sole helium line observed at high resolution with 
HST/STIS is He I A6678, which is observed in absorption for both epochs. It was hoped that we would have 
been able to determine H/He abundances from Ha and Hei, but the Hei absorption is rather insensitive to 
the precise H/He abundance ratio, and more problematic, predicted absorption strengths are significantly 
weaker than observations. 

Another important discrepancy is that the continuum fit to observations is imperfect, at both epochs; 
(a) Fen blanketing is apparently underestimated between AA1700-2300 using our determination of the iron 
content in VI (see § 3.4); (b) the models clearly fail to reproduce the observed continuum distribution either 
side of the Balmer 'jump', by up to ^20%. Hillier et al. (2000) experienced similar difficulties in their 
analysis of r\ Car (see § 3.5). 

3.3. Mass-loss properties 

In addition to deriving stellar parameters, we have simultaneously modeled the Ha profile for each epoch 
by varying the mass-loss rate, and terminal velocity. 

Using volume filling factors of ^30%, which provide a reasonable match to the electron scattering wings 
of Ha, we find an approximately constant M=4.4x 10~ 4 M Q yr _1 for both epochs, such that the stronger 
Ha emission in 1999 results from the increase in stellar temperature, rather than a change in mass-loss rate. 
Wind velocities are also very similar, ^=325 kms -1 and Uoo=290 kms -1 for 1997 November and 1999 July, 
respectively. Differences in Ha absorption can be matched with different turbulent velocities, Wturb=10 and 
40 kms -1 . Comparisons between synthetic spectra and observations are shown in Figure 6. 

An important caveat which should be kept in mind is that we have relied solely on Ha as a mass-loss 
indicator, since all other Balmer lines are unresolved. Nevertheless, low dispersion HST/STIS observations 
reveal that the increase in strength of Ha from 1997 to 1999 is not repeated for H/3 and higher Balmer 
members (Table 2). Were we instead to rely on the H/3 flux, our modeling would reveal a ~30% decrease 
in mass-loss, rather than a constant rate, between the two epochs. Such differences in behavior amongst 
members of the Balmer series are not repeated amongst other LBV-type supergiants (e.g. HDE 316285: 
Hillier et al. 1998). 

One further complication is that since we have used solely Ha to constrain the mass-loss rate, this 
is also dependent on the atmospheric H/He composition. For a solar composition the clumped mass-loss 
rates would decrease to ~ 2.5 x 10~ 4 M Q yr" 1 . Conversely, if VI is more highly enriched (e.g. H/He=l by 
number), we would need to increase the mass-loss rate estimate by a factor of ^2 to 10~ 3 M Q yr _1 . 

3.4. Metal Abundances 

Unfortunately the low spectral resolution of our observations together with the narrow, weak emission 
from most metal ions means that the determination of elemental abundances is not possible. However, the 
presence of numerous optical Fen transitions, together with a forest of metal lines in the near-UV permit 
an estimate of the iron abundance to be determined. Neglecting the effect of charge exchange would require 
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systematically higher abundances (at least by a factor of two) than we obtain here. Ideally we would like 
to have used a similar star at known metallicity to check our modeling procedure against, but no such stars 
exist with parameters close enough to VI. 

In Figure 7 we compare de-reddened ultraviolet observations of VI from 1997 November with Fell 
blanketed models at 0.20Z Q , 0.5Z Q and 1.0Z & degraded to the resolution of the STIS spectroscopy. Each 
model provides a fair match to the observed continuum depression around 2500A including the discontinuity 
at A2630. The 0.20Zq case compares very closely with observations in the AA2500-2900 region. The 0.10Z Q 
model, not shown in this Figure, is marginally worse than the 0.20^0 case. Discrepancies remain for each 
case, especially at AA1750-2300. We attribute such discrepancies to the shear complexity of the observed 
UV spectrum (Hillier et al. 2000 obtained similar problems for models of r\ Car). For example, in the case 
of the I.QZq model, 2878 lines are predicted to have equivalent widths of >0.lA in the 1200-3100A region. 
These are dominated by 2351 lines of Fen and 416 lines of Fein. In the light of model discrepancies between 
A1750-2300, a comparison with observations marginally favors the low metallicity cases, but is unable to 
provide a definitive abundance determination. 

Figure 8 compares rectified observations in the 3000-5500A region, with synthetic spectra covering the 
same range in metallicity. Aside from the Balmer series, the majority of features are again due to Fell, 
the strengths of which are very sensitive to metal content. As for the UV, an SMC-likc iron abundance 
for VI is favored by the weakness of iron lines in the AA3100-3500, AA4450-4600 and AA4900-5300 regions. 
Unfortunately, prominent P Cygni Nai emission at A5890-6 is not predicted sufficiently strong in any case, 
while comparisons in the far-red are prevented by significant CCD fringing. 

Note, however, that caution is necessary. Although we include over 700 energy levels of Fe n in our 
calculations, this ion is extremely complex, with potentially imprecise atomic data. In addition, although 
Fell is certainly the most prominent heavy element, ions of other less abundant elements will make significant 
contributions to the line blanketing at particular wavelengths (e.g. Mil). Work is currently underway to 
include such elements in our calculations. 

Further, we have based our iron abundance estimate on a comparison between models selected to match 
the 1997 November HST/STIS dataset. Were we instead to have used the 1999 July observation, we would 
have revealed a somewhat lower iron abundance, since the emission strength of optical Fe n lines (e.g. A5015) 
is weaker at that epoch. 

Nevertheless, we conclude that the heavy metal content of NGC 2363-V1 is ^0.2Z Q , in reasonable 
agreement with the previously derived nebular oxygen abundance of 12+log(0/H)=7.9, a factor of six below 
the solar value. Recall that iron is produced by supernovae, whose progenitors have lifetimes of typically 
lGyr, while oxygen is created exclusively by short-lived, high mass stars. Our study represents the most 
distant star whose iron abundance has been derived, and compared with a-elements. 

3.5. Balmer jump evidence for asymmetry? 

The shape of the Balmer jump is notoriously difficult to reproduce in extreme OB supergiants. The 
case for VI is no exception - from Figure 5 the model does not fit the continuum well on both sides of the 
Balmer jump. In fact, the model predicts a lower continuum level shortward of the Balmer edge while the 
observations clearly show an increasing flux below 3700 A. This is best illustrated for VI in Figure 9, and 
includes HST/FOS spectra of two B stars in the LMC cluster Breysacher 73 (Walborn et al. 1999). The 
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spectral shape of VI around the Balmer edge is quite similar to that of Brey 73-1C, an 09.5-Blpe star, 
classified as such because an emission component is clearly seen in the center of the Balmer absorption lines. 
The unusual flux increase from 3750 to 3650 A which is seen in both stars, has been observed in a number 
of Be stars (Kaiser 1987) and has been interpreted as evidence for the existence of a disk-like extended 
envelope of low density around the stars. The mass-loss rate of VI is much higher than that of Be stars 
and the two-component model may not be applicable in this case, but the peculiar shape of its continuum 
around the Balmer edge may indicate that its envelope is not spherically symmetric. Indeed, most nebulae 
around LBVs exhibit axisymmetric morphologies that could have been shaped by non-spherical mass loss 
(Nota et al. 1995). Spectropolarimetry of the Ha line and its adjacent continuum could be used to measure 
the degree of asymmetry, if any, of VI 's envelope. 

4. DISCUSSION 

Our results place tight constraints on the current stellar properties of NGC2363-V1. Major LBV 
eruptions are thought to coincide with an increase in the bolomctric luminosity of the star and probably 
also its mass-loss rate. We now consider whether VI has undergone a genuine eruption, and subsequently 
compare its properties with other LBVs, notably HD 5980 during its brief eruption in 1994 (Barba & Niemela 
1994). 

4.1. Is NGC 2363- VI undergoing a genuine eruption? 

Since the presence of NGC 2363-V1 was unknown prior to its recent brightening, no spectroscopic or 
photometric information is available to us in order to readily constrain its pre-outburst properties. Conse- 
quently, we consider three possibilities for its spectral type prior to 1993, namely an early O, late O, or early 
B supergiant. 

For 03 and 08 supergiant cases, we adopt a temperature and bolometric correction scale from Vacca 
et al. (1996), together with mass-loss rate estimates following the wind momentum-luminosity relation of 
Puis et al. (1996) for low metallicity (mixture of LMC and SMC) stars. For the early B supergiant case, we 
consider the possibility that VI was a dormant, hot LBV, equivalent to AG Car or R127 during their WN11 
phase (e.g. Smith et al. 1994), whose wind properties are estimated from Crowther (1997). Properties 
derived in this way are presented in Table 5. 

We first consider the possibility that VI was a hot, dormant LBV immediately prior to 1993. In this 
case, the low bolometric correction would indicate that the luminosity of VI would have been amongst the 
lowest of all known LBVs with log(L/L ) < 5.5, comparable with R110 (Stahl et al. 1990). In contrast, 
those LBVs which have been know to experience a very hot WN11 phase have luminosities in excess of 
log(L/L Q ) > 6. Nevertheless, for such a scenario, the increase in bolometric luminosity corresponds to a 
factor of approximately ten, with a mass-loss rate probably increasing by several hundred fold. We neglect 
the possibility that VI was a cooler, P Cygni-type dormant LBV since its bolometric luminosity would then 
have been unrealistically low. 

Alternatively, VI may have advanced immediately from an O supergiant to an erupting LBV. The age 
of NGC2363-B, from which VI is thought to have originated is 3-5Myr (Drissen et al. 2000), so an 08 
supergiant progenitor represents a realistic possibility. In such a case, the increase in luminosity by VI is 
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less severe (see Table 5), albeit a factor of four (or greater), while the increase in mass- loss rate is perhaps 
one thousand fold if the progenitor has wind characteristics comparable to Magellanic Cloud supergiants. 
The absence of cool dust around VI favors an O supcrgiant progenitor rather than a dormant LBV. 

The final possibility, is that VI advanced directly from an 03 supergiant. These very young (<2- 
2.5 Myr) stars have the highest bolomctric corrections of any hot massive star, prior to the Wolf-Rayet 
phase. We consider this possibility highly unlikely, given their youth relative to NGC2363-B. Nevertheless, 
for completeness, the change in stellar and wind properties would be the least severe, with potentially only 
a modest increase in luminosity, together with a factor of one hundred increase in mass-loss rate. 

The rapid brightening of VI also indicates a corresponding increase in the physical size of the star. 
Based upon a late O supergiant progenitor, the surface of VI expanded from a radius of ~ 20 to 42Oi? 
in a little over two years (1992 October to 1995 February), corresponding to an average rate of +4 kms -1 ! 
This calculation assumes that the stellar temperature of VI in early 1995 was comparable to that of 1997 
November, which appears realistic considering the similarity in photometry (Table 1). Subsequently, between 
1997 November and 1999 July, the radius of VI decreased from 420 to 315i? Q , at an average rate of —1.4 
kms -1 . These changes support some form of pulsational or dynamical instability as the cause of the present 
eruption (e.g. Guzik et al. 1999; Stothers 1999a). Calculations in our group arc underway to assess such 
possibilities. 

In summary, it appears that VI is a moderately massive star which underwent an dramatic increase 
in luminosity, radius and mass-loss rate between late 1992 and early 1995, with a probable late O-type 
supergiant or perhaps hot, dormant LBV progenitor. Estimating an initial mass for VI is extremely difficult 
since this would rely on a comparison between the luminosity of the star prior to eruption with low metallicity 
evolutionary tracks (Meynet et al. 1994). Nevertheless, an age of 4-5 Myr suggests an initial 4O-6OM star, 
with a present luminosity of log(L/L©) = 5.8 — 6.0, consistent with a late-type O supergiant precursor. 

Consequently, we suggest that the initial and current mass of NGC 2363-V1 is no more extreme than 
many other LBVs, such as AG Car. In contrast, the initial mass of the Pistol star appears to be remarkably 
high (>200M Q , according to Figer et al. 1998). 

4.2. Comparison between VI and other LBVs 

We have shown that VI is currently undergoing a giant eruption, and that its stellar and wind properties 
are exceptional. What distinguishes VI from 'dormant' LBVs, such as AG Car? How do these compare with 
HD 5980 which briefly underwent an eruption in a similar low metallicity environment during 1994? 

The well known Galactic LBV AG Car is currently undergoing huge variations in temperature (from 
8, 000-25, 000K) over timescales of decades. Its stellar properties have been studied by Leitherer et al. 
(1994) who found a fairly (homogeneous) mass-loss rate of ^10 _5 M Q yr _1 at a constant luminosity of log 
(L/Lq)=6.1. Consequently, VI is currently little more than a factor of two times more luminous than AG 
Car, which is surprisingly small. In contrast, the mass-loss rate of VI is perhaps 50 times higher, and appears 
to be the principal defining characteristic (see below). 

Are the properties of VI during its present eruption typical of other massive stars? This question is 
particularly relevant to stellar evolution models of massive stars in which, brief episodes of very intensive 
mass-loss are predicted or assumed. To address this question, we have investigated the stellar properties of 
HD 5980 during its brief eruption (Barba & Niemela 1994) based on the optical spectroscopy obtained by 
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Heydari-Malayeri et al. (1997) during the visual maximum in 1994 September when V=8.6 mag (Albert 
Jones, private communication). Based on an interstellar reddening of E(B-V)=0.07 mag towards HD 5980 
(Crowther 2000), together with an SMC distance modulus of 18.9 mag (Westerlund 1997), the absolute visual 
magnitude of the erupting star was My = —10.5 mag (all other components play a negligible contribution 
to the continuum flux at this epoch). This is remarkably similar to the current absolute visual magnitude of 
VI. However, the spectral appearance of the HD5980 eruptor was considerably different from VI at visual 
maximum. Heydari-Malayeri et al. (1997) revealed a much hotter, WNll-type appearance (Fig. 10) so it 
was intrinsically much more luminous than VI. 

We have analysed the optical spectrum of HD 5980 from Heydari-Malayeri et al. (1997), using the Hillicr 
& Miller (1998) code. Based on the usual Hci A4471, Hen A4686 and H/3 diagnostics, our derived parameters 
are listed in Table 4, revealing that the luminosity of the HD 5980 eruptor exceeded both VI and r\ Car by a 
factor 3-6, respectively. Indeed the only other LBV whose luminosity may compare favourably with HD 5980 
is the Pistol star (Figer et al. 1998), although analysis of it is severely hampered by the high reddening plus 
lack of definitive temperature constraints. A comparison between our new results for HD 5980 with those of 
Koenigsberger et al. (1998) based on observations from 1994 December (when it had visually faded by 1.6 
mag) suggests that its bolometric luminosity had decreased by a factor of five over three months! Note also 
the discrepant H/He abundance pattern resulting from analysis of the two datasets in Table 6. 

The mass-loss rate for HD 5980 during its brief eruption compares very closely with VI. It appears that 
'giant eruptions' of LBV's in low mctallicity environments exhibit representative mass-loss rates approaching 
^10 _3 M Q yr^ 1 . These differ greatly from the sole source of mass-loss rate produced in giant LBV eruptions 
known to date, namely the famous Homunculus produced by r\ Car between 1837-1960. During this eruption, 
<~ 2 — 3M Q were ejected (Davidson & Humphreys 1997) corresponding to ^10 _1 M Q yr _1 . Thus far, r\ Car 
appears to be unique in exhibiting such high rates of mass-loss, namely during its 1840's eruption. Indeed, 
although not presently in outburst, r/ Car remains a very extreme star, as illustrated in Figure 10, with a 
present luminosity and mass- loss rate somewhat higher than VI (see also Table 6). Hence, the present VI 
outburst should not be considered an equivalent to the great rj Car 1837-1860 eruption. 

5. CONCLUSIONS 

From our quantitative analysis of an HST/STIS spectrogram of NGC 2363- VI, and a photometric follow- 
up of this star since the beginning of its present eruption, we conclude as follows. VI is a Luminous Blue 
Variable currently experiencing a major eruption. Its immediate progenitor star was probably an initial 
«60M Q late O supergiant. 

Its physical parameters are exceptional: log (L/Lq) = 6.35, T* = llkK, = 420R Q for 1997 November; 
and log (L/Lq) = 6.4, T* = 13kK, i?* = 315i? Q for 1999 July. Therefore, the bolometric luminosity of VI is 
exceeded by only, in ascending order, r\ Car (log L/Lq — 6.7, Hillier et al. 2000), the Pistol star (log L/Lq 
= 6.9±0.3, Figer et al. 1998) and HD 5980 during eruption (log L/Lq = 7.1). 

Analysis of high resolution Ha observations indicate that the wind is moderately clumped, with M ~ 
4.5xlO~ 4 M Q yr _1 and Voo ~300 kms -1 . During its brief eruption in 1994, HD 5980 possessed a similar 
mass-loss rate. Therefore, the mass-loss rate of ~ 1O _1 M yr _1 exhibited by r\ Car during its giant eruption 
in 1840's (Humphreys & Davidson 1994) is thus far unique. 

The presence of many iron features in the spectrum permits an estimate of its heavy metal content. 
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Although discrepancies remain, especially at UV wavelengths, our study reveals a SMC-like Fe content 
for NGC2363-V1, which is in good agreement with previously derived oxygen abundances for NGC2363. 
Consequently, the heavy metal enrichment of NGC 2363, produced by Type I Supernovae, whose progenitors 
are believed to have a lifetime of typically 1 Gyr, appears to mimic that of light elements, such as oxygen, 
which are created exclusively by short-lived (lOMyr) massive stars. 

Observationally, spectroscopic and photometric monitoring of VI and other LBVs will continue, while 
theoretical studies are underway to investigate the cause of the outbursts experiences by VI and HD 5980, 
whether pulsationally, dynamically or rotationally driven. 

We would like to thank Mohammad Heydari-Malayeri for providing the outburst spectrum of HD 5980, 
Albert Jones for providing his photomety of HD 5980, and Kris Davidson and his HST/STIS observing team 
for providing the spectrum of Eta Carina in Figure 10. This investigation was funded in part by the Natural 
Sciences and Engineering Research Council of Canada, by the Fonds FCAR of the Government of Quebec 
and by Universite Laval (LD, JRR and CR). PAC acknowledges financial support from the Royal Society. 
We are especially grateful to the staff of the now defunct Royal Greenwich Observatory for obtaining service 
imaging of NGC 2363. The William Hcrschel Telescope is operated on the Island of La Palma by the Isaac 
Newton Group in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de 
Canarias. 
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Table 1. Photometric variability of NGC 2363-V1 since 1996 



Epoch 


Observatory 


F170W 


U 


B 


V 


R 


F1042M 


1996 Jan 


HST/WFPC2 






17.82 


17.88 






1997 Nov 


HST/STIS 


17.25 


16.79 


17.55 


17.43 


17.16 




1998 Jan 


WHT/aux 




16.13 


17.53 


17.66 


17.33 




1998 Mar 


HST/WFPC2 


17.23 


16.72 




17.52 




17.15 


1998 Dec 


HST/WFPC2 


17.09 


16.79 




17.55 




17.20 


1999 Jul 


HST/STIS 


16.83 


16.85 


17.72 


17.58 


17.20 




1999 Nov 


HST/WFPC2 


16.89 


16.89 




17.80 




17.56 


2000 Feb 


HST/WFPC2 


16.98 


16.72 




17.76 




17.50 



Table 2. Equivalent widths of the Balmer lines in VI 



Line 


W x (1997) 


W x (1999) 


Ha 


250 


300 


H/3 


45 


35 


H 7 


15 


11 


RS 


6 


5 



Note. — The equivalent width, in 
units of A, is for the emission com- 
ponent only. 
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Table 3. Equivalent widths of the the strongest UV lines in VI and NGC 2363 



Line 


W x (VI) 

(A) 


W x (knot A) 


W x (knot B) 




(A) 


(A) 


Cm A1175 


2.4 






Sin A1260 


2.8 


0.9 


1.6 


Sin A1265 


1.8 






Oi/Sin A1303 


4.1 


2.2 


2.3 


Sin A1309 


1.0 






Cn A1334 


2.5 


1.6 


1.6 


Si iv A1393 


2.0 


0.6 


0.6 


Si iv A1402 


2.3 


0.5 


1.1 


Sin A1527 


2.1 


0.8 


0.8 


Sin A1533 


1.9 


<0.2 


<0.1 


Civ A1549,51 


2.8 




2.1 



- 18 - 



Tabic 4. Summary of model atoms used in radiative transfer calculations 



Ion 


N s 


N F 


-^Trails 


Details 


Hi 


10 


30 


435 


n <30. 


He i 


27 


39 


315 


n <14 


He ii 


5 


5 


10 


n <5 


Ni 


44 


104 


855 


nl <5f2F° 


Nil 


23 


41 


144 


nl <2p3d 1 P° 


Nm 


8 


8 


11 


nl <2p 32 P° 


Nai 


18 


53 


507 


n <12 


Mgll 


18 


45 


362 


n <10 


Alii 


26 


44 


171 


nl <5d!D 


Aim 


12 


20 


67 


n <5 


Can 


17 


58 


463 


n <11 


Sin 


31 


49 


281 


nl <3s 2 9h2H° 


Sim 


12 


20 


42 


nl <3s4p 1 P° 


Fen 


134 


827 


21907 


nl <3d 5 ( 6 S)4p 2 ( 3 P) 4 P 


Fein 


69 


607 


6620 


nl <3d 5 ( 4 D)6s 3 D 



Note. — Including full levels (Np), super levels (Ns) 
and total number of transitions (-/Vxrans), following Hillier 
& Miller (1998). 
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Table 5. Estimates of stellar properties for NGC 2363-V1 prior to outburst 



Sp. Type 




i?* 


log L* 


logM 






kK 


Rq 


L Q 


Mqyt- 1 


kms -1 


031 


50.7 


18 


<6.27 


-5.3 


3200 


081 


35.7 


22 


<5.85 


-6.2 


2000 


B0.5I 


25.6 


27 


<5.46 


-5.8 


250 



Note. — We use here the temperature and bolometric 
correction scale from Vacca et al. (1996), plus the wind 
momentum-luminosity relation of Puis et al. (1996) for 
low metallicity (LMC and SMC) OB stars, except for the 
early B supergiant for which typical LBV wind properties 
are shown (Crowther 1997). 



Table 6. Stellar parameters of VI and of other well known LBVs 



Star 


Epoch 


T/kK 


R/Rq 


log L» 


logM/V? 




H/He 


My 


Reference 






TR,o SS = 10 


TRoss = 2/3 


L 


Meyr- 1 


kms -1 




mag 




NGC 2363-V1 


Nov 1997 


11 


590 


6.35 


-3.1 


325 


4.0? 


-10.4 


This work 




Jul 1999 


13 


675 


6.4 


-3.1 


290 


4.0? 


-10.3 


This work 


HD 5980 


Sep 1994 


23 


280 


7.05 


-3.1 


500 


0.7 


-10.5 


This work 




Dec 1994 


35 


130: 


6.5 


-3.0 


600 


2.5 


-8.6 


(1) 


rj Car 


Mar 1998 


27 


881 


6.7 


-2.5 


500 


5? 


_? 


(2) 



Note. - - We estimate that a volume filling factor of ,/=0.3 is appropriate for VI, with /=0.1 for 77 Car and 
HD5980. References: (1) Koenigsberger et al. (1998); (2) Hillier et al. (2000). 
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Fig. 1. — V-band light curve of NGC 2363-V1 over the past 7 years. The absolute magnitude scale on the 
right assumes a distance of 3.44 Mpc and E(B — V) = 0.06 (see text). 
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Fig. 2.— Flux calibrated HST/STIS spectrogram of VI from 1999 July (solid) and 1997 November (dotted). 
The ripples longward of 7500 A are due to fringing on the CCD. 
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Fig. 3. — High resolution HST/STIS spectrogram of VI (grating G750M), revealing a P Cygni Ha profile. 
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Fig. 4.— Rectified spectra of NGC 2363-V1 in 1999, NGC 2363- A, NGC 2363-B and a representative Galactic 
B5.5I spectrum (see text). 
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Fig. 5. — Comparison between de-reddened HST/STIS datasets (solid) from 1997 November (upper panels) 
or 1999 July (lower panels) with non-LTE model predictions (dotted). 



- 25 - 




Fig. 6. — Comparison between rectified Ha profiles (solid) and synthetic spectra (dotted). 
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Fig. 7.— Comparison between de-reddened STIS/MAMA datasets of VI with line blanketed UV energy 
distributions, for O.2OZ (lower), 0.5Z Q (middle), I.QZq (upper), for clarity successively offset by 2.5xl0~ 15 
erg cm -2 s _1 A. 
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Fig. 8. — Comparison between rectified, optical STIS datasets of VI with line blanketed energy distributions, 
for 0.20Z Q (lower), 0.5Z Q (middle), I.OZq (upper), for clarity successively offset by a continuum unit. 
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Fig. 9. — Comparison between VI and two LMC B stars in the Breysacher 73 cluster, 1C (09.5-Blpc) 
and # 6 (B0.2V), in the vicinity of the Balmer jump. 
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Fig. 10. — Comparison between rectified, optical datasets of VI with HD5980 during outburst in 1994 
September, plus rj Car in 1998 March (HST/STIS data from Davidson et al. , in preparation). 



